MULTIVARIATE ANALYSIS OF THE VEGETATION IN A TWO-DESERT INTERFACE 


A. A. El-Ghonemy 1 , A. Wallace-, and E. M. Homney : 


Abstract.— This report further describes the distribution and ecological characteristics of the natural vegetation 
at the Mojave Desert-Great Basin Desert interface. The region studied is one of extraordinary biological interest 
because of its geographic location straddling the boundaries of two large deserts of the western United States, and 
because of the kind and manner of its past land use (atmospheric and underground testing of nuclear devices). The 
present analysis determines the magnitude of variations in the phytosociological structure in this region and eval¬ 
uates some relationships between its vegetation and environment. Vegetation and soils were sampled in 66 stands 
representing many possible physiographic variations. Relative density and relative coverage were determined for 
each perennial species and summed to provide an estimate of its importance value (I.V.). Importance values were 
used to ordinate stands to provide a synthesis of the phytosociological data and to portray the compositional relation¬ 
ships of species. The results of this study indicate that the area is dominated bv several interrelated vegetational 
groupings. Correlations between the vegetational groups and the different environmental variables indicate that the 
distributional pattern of the vegetation is controlled largely by soil physical properties, salinity, and fertility levels. 


The landscape of the Nevada Test Site is 
one of the most intensively studied and best 
understood deserts in the United States. Its 
potential is unique for studies critical to bet¬ 
ter understanding of arid lands. Previous phy¬ 
tosociological studies in this area are largely 
descriptive (Wallace and Romney 1972, 
Romney et al. 1973, Beatley 1976). 

The objective of this study was to use some 
multivariate methods to analyze sociological 
relations among plant communities of the 
natural vegetation on the Nevada Test Site. 
This study is closely related to that pre¬ 
viously carried out by the authors (El-Gho¬ 
nemy et al. 1980), in which an account is giv¬ 
en on the location, physiography, climate, 
vegetational groupings, and community di¬ 
versity. 

Materials and Methods 

Procedural details involving selection of 
stands and sampling techniques for soil and 
plants at 66 sites (Table 1) have been report¬ 
ed by Wallace and Romney (1972) and Rom¬ 
ney et al. (1973). For treatment of data, one 
classification and two ordination techniques 
were used to analyze a data matrix consisting 
of the importance values for each of the pe¬ 


rennial species encountered in each of the 
stands. 

The classification technique involves the 
unweighted pair-group agglomerative clus¬ 
tering, using arithmetic averages to compute 
the similarity between a cluster and a stand 
which is a candidate for entry into a cluster 
(Sneath and Sokal 1973). The Euclidean dis¬ 
tances (Ed) were used as the measure of sim¬ 
ilarity between stands. 

The first ordination technique is that of 
Wisconsin (Gray and Curtis 1957) as modi¬ 
fied by Beals (1969). The raw data were nor¬ 
malized by row and column, and interstand 
similarities were calculated using the formula 
2w -s- (a + b), where w represents the sum 
of the smaller values for common species; a 
and b represent the sum of all species in 
stands A and B, respectively. The maximum 
dissimilarity value was set equal to the max¬ 
imum similarity value found in the similarity 
matrix. Ten percent of the i’th axis was 
searched to locate the second end stand for 
the (i + l)’th axis. 

The second ordination technique involves 
principal component analysis of the matrix of 
interstand correlation coefficients (Sneath 
and Sokal 1973). Eigenvectors (normalized to 
the eigenvalues) were not rotated. 
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Results 

1. Classification of the Vegetation Data 

(a) The clustering units.— The program 
cluster was used and Ed was selected as a 
measure of similarity. It should be mentioned 
that with this method of classification the 
stands are clustered into cells regardless of 
whether they form discrete groups in nature 
or whether they are merely parts of a contin¬ 
uum. The dendrogram shown in Figure 1 is 
derived from this cluster analysis. Because 
the paired and grouped stand clusters in the 
dendrogram result in linkages at various lev¬ 
els of similarity, an ecologically meaningful 
classification is not automatically indicated. 
Classification can be obtained, however, by 
setting more or less arbitrary threshold val¬ 
ues. Threshold values used were set at 10, 15, 
20, and 29 Ed. These are indicated by hori¬ 
zontal dash lines on the dendrogram. 

At Threshold Line 4 there are three main 
clusters. Cluster 1 links together 58 stands (1 


through 62) at Ed distance of about 28. Clus¬ 
ter II links together 5 stands (21 through 24) 
at Ed of 20. Cluster III links together stands 
33, 39, and 35 at Ed of about 16. 

At Threshold Line 3 clusters II and III re¬ 
mained unaltered, but cluster I became dis¬ 
tinguishable as three subclusters: subcluster 
IA that links together 50 stands (1 through 
65), subcluster IB that links together stands 
54 and 62. At this level of similarity stands 
19, 31, and 11 became so dissimilar to other 
stands that they remain rather isolated. 

At Threshold Line 2 the subcluster IA be¬ 
came distinguishable as five vegetational 
groupings. The first grouping, IAa, links to¬ 
gether 11 stands (3 through 34). The second 
grouping, IAb, links together 12 stands (6 
through 51) at Ed of 12. The third grouping, 
IAc, links 8 stands (4 through 15) at Ed of 
about 14. The fourth grouping, lAd, links 6 
stands (12 through 48) at Ed of about 14.5. 
The fifth grouping, IAe, links together stands 
50, 64, 66, and 65 at Ed of about 16. At this 
level of similarity many individual stands or 



Fig. 1. Dendrogram resulting from the application of the agglomerative clustering analysis. The dotted lines de¬ 
note the levels at which the dendrogram yields different vegetational groupings. 
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Tshle 1. Relative dominance of major species in each of the 66 stands (those 2 percent or more lor a iyven stand 
(See Table 2 for plant abbreviations.! 
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Table 1 continued. 
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couples of stands become mostly dissimilar to 
the well-defined clusters so that they remain 
isolated (Fig. 1). 

Classification at Threshold Line 1 resulted 
in many fragmentary units of limited general- 
izable value. 

Sociological significance of the vegetation 
groupings.— For purposes of discussion each 
of the clusters identified, irrespective of its 
hierarchal level on the dendrogram, was 
called a vegetational grouping and named af¬ 
ter the most abundant species, that is, the 
species with the highest average importance 
value. 

Table 2 includes the average importance 
values for the different species in the various 
vegetational groupings. Inspection of this 
table gives the following explanation of the 
results of the cluster analysis. 

Vegetational Grouping lAa—Larrea triden- 
tata (Sesse & Moc. ex DC) Cov.: This group¬ 
ing is represented by stands from Frenchman 
Flat (7 stands), Mercury Valley (2), Jackass 
Flats (1), and Rock Valley (1). The clustering 
seems to be based on the presence of L. tri- 
dentata as a leading dominant in all the 
stands in the cluster. The average importance 
value of L. tridentata (Table 2) is 93.4 (out of 
200). The associated species are generally of 
minor importance, except Ambrosia dumosa 
(A. Gray) Payne (I.V. = 23.8). Stand 36 (Fig. 
1), although dominated by L. tridentata (I.V. 
= 112), is the last to join the cluster. This is 
primarily due to the presence of Atriplex ca~ 
nescens (Pursh) Nutt, in substantial amounts 
(I.V. = 55). 

Vegetational Grouping IAb—A. dumosa: 
This cluster does not appear to be a very nat¬ 
ural unit, and on the basis of the subleading 
dominant species there would be grounds for 
the recognition of three smaller clusters (Fig. 
1). The most influential species responsible 
for the segregation of this vegetational 
grouping into smaller clusters are L. triden¬ 
tata (in stands 6, 30, 37), Krameria parvifolia 
Benth and Coleogyne ramosissima Torr. (in 
stands 7, 8, 55, 56), Grat/ia spinosa (Hook.) 
Moq. (in stands 40, 44, 49, 51), and Onjzopsis 
hymenoides (Roem. & Schult.) Ricker (in 
stand 60). The average importance value for 
A. dumosa in this grouping is 77. Associated 
species of pronounced significance are L. tri¬ 


dentata (I.V. = 37) and K. parvifolia (I.V. = 
29.4). 

Vegetational Grouping 1 Ad—Transitional: 
Stands of this grouping also show little re¬ 
semblance and mostly fuse in pairs above an 
Ed of about 10. Dominance is shared by 
many species, among them Lycium atider- 
sonii A. Gray (I.V. = 22.2), Ephedra ncva- 
densis S. Wats. (I.V. = 22.1), O. hymenoides 
(I.V. = 17), and C. ramosissima (I.V. = 15). 

Vegetational Grouping IAe—G. spinosa: 
The relatively high similarity between stands 
64 and 68 (Fig. 1) is not only due to the dom¬ 
inance of G. spinosa, but also due to the pre¬ 
ponderance of L. tridentata, E. nevadensis, 
and C. ramosissima. The average importance 
value of G. spinosa in this grouping is 71, fol¬ 
lowed by 23 for L. andcrsonii and 22.5 for C. 
ramosissima. 

Vegetational Grouping IB— Lycium 
shockleyi A. Gray (L. rickardii) C. H. Mull.: 
All stands of this cluster are also dominated 
by L. shockleyi, with an average importance 
value of 88. Other important species are L. 
tridentata (I.V. = 35) and Atriplex confer- 
tifolia (Torr. & Frem.) S. Wats. (I.V. = 27). 

Vegetational Grouping IC—C. ramosis¬ 
sima: All of the five stands representative of 
this cluster are dominated bv A. canescens. 
Stands 28 and 38, which fuse at very high 
similarity levels, are overwhelmingly domi¬ 
nated by A. canescens (I.V. = 188 and 185, 
respectively). The average importance value 
for A. canescens in this cluster is 150. Species 
of some importance are Stanleya pinnata 
(Pursh) Britt. (I.V. = 18) and Lycium palli¬ 
dum Miers (I.V. = 8.4). 

Vegetational Grouping III—A. confer- 
tifolia: This is the last cluster to join the den¬ 
drogram, and it enters at a very low level of 
similarity (Ed = about 33). All stands are 
dominated by A. confertifolia. Stand 33 (I.V. 
= 200) Rises with 39 (I.V. = 193) at the 
highest recorded level of similarity (Ed = 
0.8). Stand 35, which is also dominated by A. 
canescens, has a high (45.6) importance value 
for A. canescens as well as for Ceratoides la- 
nata (Pursh) J. T. Howell (I.V. = 22). This 
results in its fusion with stands 33 and 39 at a 
relatively low similarity level (Fig. 1). Other 
species in this grouping are of little signifi¬ 
cance. 
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2. Ordination of the Vegetation Data 

(a) Ordination of Stands : According to the 
Brav and Curtis (1957) technique, the two-di¬ 
mensional ordination of stands (Fig. 2) has re¬ 
sulted in three distinct hyperspheres which 
correspond to the vegetation groupings iden¬ 
tified at Threshold Line 4 on the dendrogram 

I derived from the cluster analysis and super¬ 
imposed on the ordination plane. Smaller 
vegetational groupings identified at higher 


similarity levels (Threshold 2) are mostly in¬ 
terconnected on the ordination plane (see 
dashed lines), but are still distinguishable in 
the form of successive groups of stands segre¬ 
gated along the primary X axis. 

In the application of principal component 
analysis (PCA), five components or axes were 
extracted that account for 65.3 percent of the 
total variation (Table 3). Plotting of stand 
scores (Figs. 3 and 4) on the two axes that 
showed a greater number of significant corre- 


Table 2. Average importance values for 35 perennial species in the different vegetational groupings (see Fig. 1). 
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Larrea tridentata (Lt)°° 

93.4 

39.2 

30.3 

41.8 

13.5 

38 

24.0 

_ 


Ambrosia dumosa (Ad) 

23.8 

77.4 

15.8 

14.5 

14.0 

2 

0.3 

0.5 

1.3 

Oryzopsis hymenoides (Oh) 

7.0 

7.2 

10.1 

1.0 

3.3 

3 

_ 

3.1 

2.2 

Ceratoidcs lanata (Cl) 

5.0 

3.6 

23.7 

4.0 

2.9 

_ 

0.6 

1.0 

9.3 

S phaeralcea ambigua (Sa) 

2.6 

1.2 

2.0 

1.2 

- 

2 

3.3 

2.2 

- 

Lt/cium andcrsonii (La) 

16.0 

5.3 

17.4 

38.0 

23.0 

_ 

11.0 

2.4 


Mirabilis pudica (Mp) 

0.1 

- 

- 

2.6 

- 

1 

- 

3.0 

_ 

Hymenoclca salsola (Hs) 

0.8 

— 

2.5 

11.2 

5.1 

1 

8.1 

_ 

2.7 

Acamptopappus shockley i (As) 

6.8 

7.4 

44.9 

- 

_ 

21 

_ 

_ 

_ 

Ephedra funerea (Ef) 

2.0 

1.1 

2.3 

- 

- 

4 

- 

- 

- 

Psorothamnus fremontii (Pfr) 

5.6 

0.7 

0.1 

9.2 






Krameria parvifolia Kp) 

5.2 

29.4 

3.9 

22.2 

_ 

5 

_ 

_ 
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Yucca schidigera (Vs) 

3.3 

- 

0.8 

— 

_ 

_ 

_ 

_ 

_ 

Yucca brevifolia (Yb) 

1.8 

- 

- 

5.0 

— 

2 

_ 

_ 

_ 

Lycium pallidum (Lp) 

3.1 

3.1 

9.0 

10.3 

- 

- 

- 

8.4 

- 

Lepidium fremontii (Lf) 

0.3 

0.4 

0.5 

1.9 

_ 

_ 

_ 

_ 


Ephedra nevadensis (En) 

5.5 

6.8 

3.9 

22.1 

20.0 

_ 

6.3 

_ 

_ 

A tenodora spinescens (Ms) 

- 

1.3 

0.5 

0.5 

- 

4 

0.8 

_ 

_ 

Craijia spinosa (Gs) 

11.1 

10.3 

37.4 

17.0 

71.5 

— 

- 

0.9 

0.6 

Coleogyne ramosissima (Cr) 

- 

4.5 

- 

15.0 

22.5 

6 

137.0 

- 

- 

Machacranthera tortifolia (Mt) 

_ 

1.0 

0.6 

5.1 

_ 


2.4 



Cactus spp. (Ca) 

1.0 

- 

- 

- 

0.2 

- 

0.4 

_ 

_ 

S tipa spcciosa (Ssp) 

- 

- 

0.9 

- 

22.0 

_ 

0.8 

_ 

_ 

Tetradymia axillaris (Ta) 

- 

- 

1.0 

0.3 

0.4 

_ 

3.3 

_ 

_ 

Atriplex confertifolia (Aco) 

- 

1.0 

- 

- 

- 

23 

0.4 

12.8 

188.0 

Atriplex canescens (Aca) 

5.0 

1.0 

_ 

_ 

_ 

_ 

_ 

150.0 

15.0 

S tanleya pinnata (Sp) 

- 

- 

_ 

0.3 

_ 

_ 

_ 

18.0 

_ 

Hilaria rigida (Hr) 

— 

0.1 

0.3 

_ 

_ 

_ 

_ 

_ 

_ 

Lijcium shockleyi (Ls) 


- 

- 

- 

- 

88 

- 

_ 

_ 

Artemisia spinescens (Asp) 

- 

- 

0.6 

- 

- 

- 

- 

- 

- 

Primus fasciculata (Pf) 

_ 

_ 

_ 

_ 

_ 

_ 

_ 



Salazaria mexicana (Sm) 

_ 

_ 

_ 

_ 

_ 

1 

_ 

_ 

_ 

Thamnosma montana (Tin) 

0.2 

_ 

_ 

_ 

_ 

_ 

_ 

_ 

_ 

Haplopappus cooperi (He) 

- 

— 

0.6 

_ 

_ 

_ 

_ 

_ 

_ 

Psorothamnus polyadenitis (Pp) 

1.0 

- 

- 

- 

- 

- 

- 

- 

- 


’Transitional 

Letters in parentheses indicate abbreviations for species’ names. Table headings correspond to these abbreviations. 
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lations with environmental variables, viz.. 
Axes I and III, show that the application of 
PCA technique has resulted in much better 
segregation of stands into groups comparable 
with those derived from the clustering tech¬ 
nique at Threshold Line 2 on the dendro¬ 
gram. 

The two vegetation groupings of A. canes¬ 
cens and A. confertifolia are separated at the 


negative side of the first axis. Separation of 
the two groupings from one another is affect¬ 
ed by the third axis. Stand 35, which joins the 
A. confertifolia grouping at a relatively low 
level of similarity (Fig. I), appears on the or¬ 
dination plane to be more associated with the 
A. canescens grouping. The groupings of L. 
tridentata and A. dwnosa , which have more 
or less overlapping scores on the first axis, are 



x 

Fig. 2. Ordination of the stands on the first axes derived from the Wisconsin ordination (Bray 6c Curtis) technique. 
The classification of stands according to the cluster analysis has been superimposed on the ordination diagram. Dash¬ 
ed lines connect some of the vegetational groupings. Empty circles are stands of L. tridentata (IAa), crossed circles 
are stands of A. dwnosa (IAb), vertically half-blocked circles are stands of C. spinosa (IAe), horizontally half-blocked 
circles are stands of A. canescens (II), large blocked circles are stands of A. shockleyi IAc), blocked triangles are 
stands of C. ramosissima (IC), empty triangles are stands ol L. shockleyi (IB), crosses are stands of A. confertifolia 
( 111 ), and small blocked circles are transitional stands (OAd). 
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clearly separated on the third axis. Along this 
axis A. dumosa scores positive values, but 
those of L. tridentata are negative. Along the 
first axis the groupings of L. tridentata , A. 
shockleyi, L. shockleiji, and C. ramosissima 
are also separated from one another. The 
grouping of G. spinosa is primarily separated 
from other groupings along the third axis. It 
is also apparent (Fig. 4) that the groupings of 
L. shockleyi and C. ramosissima occupy a 
more or less central position between other 
groupings. 

(b) Behavior of species along the environ¬ 
mental gradients : The behavioral pattern of 


eight common species as expressed by their 
importance value is represented separately 
on the ordination derived from the first and 
third axes of the PCA (Fig. 5:A-H). Larrea 
tridentata and A. dumosa attain their high 
importance values at the high and medium 
positions of the positive side of the axis. 
Along the third axis the two species behave 
rather differently: the high values of A. du¬ 
mosa are on the positive side and the high 
values of L. tridentata are on the negative 
side. Grayia spinosa and L. andersonii also 
exhibit definite patterns, with their high im¬ 
portance values at the medium position on 


nr 
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Fig. 3. Plotting of stands on axes I and III from the principal component analysis (PCA). 
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the first axis. Acamptopappus shocklexji shows 
its high values at medium to lower positions 
on the positive side of the first axis. On the 
third axis most of the high values are on the 
negative side. More definite patterns are 
those of A. canescens and C. ramosissima. 
They have attained their high importance 
values at the extreme positive and the ex¬ 
treme negative end of the third axis, respec¬ 
tively. Atriplex confertifolia attained its high 
importance values at the central positions of 
the first axis and at medium positions with 
regard to the positive side of the third axis. 

(c) Ecological significance of plujtosociolo- 
gical gradients : Simple correlation 

coefficients between the phytosociological 
gradients represented by the five axes ex¬ 
tracted from the PCA and the various envi¬ 
ronmental variables are given in Table 3. 
Most of the correlations are extremely low. 
The lowest negative correlation (-0.39) is 
that between the first axis and field capacity 
for water relations in soil. This axis is also sig¬ 
nificantly correlated with both sodium (nega¬ 
tive) and iron (positive). Axis 2 shows no sig¬ 
nificant correlations with any of the variables 
studied. Axis 3 shows significant correlations 
and most of the variables correlated signifi¬ 
cantly with the first axis, but with opposite 
trends. This axis also shows significant corre¬ 
lations with both potassium (positive) and ni¬ 
trogen (negative). Axis 4 shows significant 
positive correlations with electrical con¬ 
ductivity and potassium, and axis 5 shows sig¬ 


nificant positive correlations with soil mois¬ 
ture retention capacity, electrical 
conductivity, calcium, magnesium, copper, 
and nitrogen. 

It is apparent from these correlation stud¬ 
ies that the segregation of the vegetation cov¬ 
er into distinct groupings on the ordination 
plane is largely attributed to variations in soil 
properties. The vegetational grouping of A. 
canescens and A. confertifolia occupy sites 
poor in phosphorus, organic matter, and ni¬ 
trogen, but rich in sodium, potassium, cop¬ 
per, and percent of clay. The grouping of A. 
dumosa reflects sites rich in phosphorus, iron, 
and to some extent in sodium, but poor in ni¬ 
trogen. The Grayia spinosa grouping oc¬ 
cupies sites decidedly poor in sodium, potas¬ 
sium, and fine particles, but rich in nitrogen 
and iron. The groupings of L. shockleyi and 
C. ramosissima occupy sites with more or less 
intermediate soil characteristics. 

(d) Correlation among species and species 
ordination : The spatial pattern of one species 
may be modified by another. This leads to 
the question of interspecific correlation. The 
causes of these correlations are, however, 
varied. The most common cause is, no doubt, 
the mutual response to varying environments. 
There also are interactions between species 
that do not involve independent environmen¬ 
tal factors (allelopathic effect, competition, 
or amelioration or degradation of environ¬ 
mental conditions). It is, however, difficult to 
reach firm conclusions as to the cause of cor- 


Table 3. Simple linear correlation coefficients (r) between five principal components and the various environmen¬ 
tal parameters 0 


Proportion of Accumulated 

Parameter total variance variance F. Lime 


PCA axes 

Percent 

Percent 

Cap. 

PH 

EC 

Percent 

1 

31.1 

31.1 

-0.39° 00 

-0.13 

0.04 

-0.20 

2 

11.1 

42.2 

0.05 

0.21 

0.08 

0.21 

3 

8.1 

50.3 

0.10 

0.07 

0.09 

-0.06 

4 

7.9 

58.2 

0.22 

-0.03 

0.27° 

-0.07 

5 

7.1 

65.3 

0.38° 00 

0.07 

0.31 00 

0.08 


*A single asterisk denotes a significant correlation at 5 percent probability level. Double and triple asterisks denote a significant correlation at 1 percent 
and 0.1 percent probability levels. 

F. Cap. is water retention at field capacity. 

EC is electrical conductivity. 

Cations are exchangeable. 

Fe, Zn, and Cu are from 0.005 M DTPA extract. 

N is total soil nitrogen. 
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Fig. 4. Ordination of stands on axes I and III from the principal component analysis. The classification of stands 
according to the agglomerative clustering analysis has been superimposed on the ordination diagram. For identi¬ 
fication of groupings see Fig. 2. 


Table 3 continued. 


Na 

K 

Ca 

plus 

Mg 

Na 

P 

Mg/g 

Fe 

Zn 

Cu 

N 

Elevation 


me/100 g 


% 

H g/g DTPA ext. 

% 

-0.34° ° ° 

-0.19 

-0.17 

-0.30° * 

-0.32* *° 

0.31°° 

-0.06 

-0.28° 

0.10 

0.18 

0.12 

0.08 

-0.10 

0.15 

-0.17 

-0.21 

0.05 

-0.22 

-0.13 

0.03 

0.35** 

0.28° 

-0.14 

0.34° ** 

-0.16 

_o 37 o.o 

-0.21 

-0.15 

-0.31 ° 0 

-0.10 

0.23 

0.34° °° 

0.17 

0.17 

0.05 

0.07 

0.14 

-0.21 

-0.10 

-0.1 

-0.12 

0.21 

0.25° 

-0.20 

0.18 

0.11 

0.13 

0.31 00 

0.27° 

0.03 
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relations by simply observing the spatial dis¬ 
tribution of the two species in nature; how¬ 
ever, if neither species separately shows any 
patterning but the two random distributions 
are coincidental (or countercoincidental), a 
direct relationship between the species seems 
the most likely explanation (Goodall 1970). 

In Table 4 a partial simple linear correla¬ 
tion matrix is given for 35 common species 
showing the positive and negative relation¬ 


ships present. The species constellation based 
on correlation values is illustrated in the form 
of a three-dimensional diagram in Figure 6. 
The components involved in the construction 
of this diagram were extracted using the 
principal component analysis of the matrix of 
interspecific coefficients (Sneath and Sokal 
1973). Five groupings of species are appar¬ 
ent—a central group and four peripherals. 
The arrangement of species in this diagram 
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Fig. 5. (A-H). The behavior of eight species as expressed by their importance values along the phvtosociological 
gradients expressed by the first and third axes of the principal component analysis. The importance value for each 
are superimposed upon a common ordination plane. High, medium, and low values are indicated where appropriate. 
A is L. tridentata , B is A. dumosa , C is G. spinosa, D is L. andersonii, E is A. shockleyi, F is A. canescens, G is C. 
ramosissima, and H is A. confertifolia. 
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reflects positive correlations in the main, the 
general position of each grouping of species 
being related to the negative correlations also 
present. The following are the five groups 
identified: 

Group 1 (upper left-hand side of the dia¬ 
gram): Graijia spinosa, L. andersonii, L. palli¬ 
dum, Stipa speciosa. 

Group 2 (lower left-hand side): Ambrosia 
dumosa , Hilaria rigida (Thurb.) Benth. ex 
Scribn., Yucca schidigcra Roezl ex Ortgies, 
Ephedra funerea Cov. & Mort., Salazaria 
mexicana Torr., E. nevadensis, K. parvifolia. 


Fig. 5 continued. 



nr 



Group 3 (upper right-hand side): Hijme - 
noclea salsola Torr. & Gray, Tetradymia ax¬ 
illaris A. Nels., A. shockleyi, Machaeranthera 
toiiifolia A. Gray, Sphaeralcea ambigua A. 
Gray. 

Group 4 (lower right-hand side): Yucca 
brevifolia Engelm. in Wats., Thamnosma 
montana Torr. & Frem., Psorothamnus fre- 
montii (Torr.) Barneby, Primus fasciculate 
(Torr.) A. Gray, Cactus species. 

Group 5 (central): Psorothamnus poly¬ 
adenitis (Torr.) Rydb., A. canescens, A. con- 
fertifolia, Mirabilis pudica Barneby, O. hiyme- 



I 
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noides , Mendora spinescens A. Gray, L. 
tridentata, S. pinnato, P. fremontii, Hap- 
lopappus cooperi (A. Gray) Hall, L. shockleyi, 
C. lanata, Artemisia spinescens D. C. Eat., C. 
ramosissima. 

It is equally clear that none of these five 
groups is an isolated entity; each is linked to 
the adjacent group by correlations of differ¬ 
ent magnitude between member species of 
the representative groups or through inter¬ 
mediate species. 

The deduction that these correlated groups 
of species represent significant communities 
in nature may not necessarily hold. The very 
obvious group on the lower right-hand side of 
the diagram (Group 4) includes species of mi¬ 


nor significance in vegetation structure (very 
low importance values), and consequently 
such a group is not detectable in nature. On 
the other hand, group 1, dominated bv G. 
spinosa, L. andersonii, S. speciosa, and C. 
ramosissima , is well defined in nature and 
floristically structurally comparable with the 
G. spinosa grouping previously defined by 
stand classification and ordination tech¬ 
niques. 

Discussion 

The data obtained by measuring the 
amount of a species in each of n stands can 
be represented by a scatter diagram of n 


Table 4. Fart of simple linear correlation matrix for 35 perennial specimens for the northern Mojave Desert show¬ 
ing positive and negative species’ relationships. 


Acamptopappus shockleyi (As) 
Atriplcx confertifolia (Aco) 
Ephedra funerea (Ef) 

Ccratoidcs lanata (Cl) 

Ambrosia dmnosa (Ad) 
Krameria parvifolia (Kp) 

Larrea tridentata (Lt) 

Lycium andersonii (La) 
Menodora spinescens (Ms) 
Yucca schidigcra (Vs) 

S phaeralcea ambigua (Sa) 
Ephedra ncvadensis (En) 
Grayia spinosa (Gs) 

Lepidium fremontii (Lf) 
Oryzopsis hymenoides (Oh) 
Machaeranthera tortifolia (Mt) 
Fsorothammis fremontii (Pfr) 

Coleogync ramosissima (Cr) 
Lycium pallidum (Lp) 
llilaria rigida (Hr) 
llymenoclea salsola (Us) 

S tanlcya pinnato (Sp) 

Stipa speciosa (Ssp) 

Tetradymia axillaris (Ta) 
Mirabi/is pudica (Mp) 

Atriplcx canesccns (Aea) 
Artemisia spinescens (Asp) 
Fsorothammis polyadenitis (Pp) 
Cactus sp. (Ca) 

Frunus fasciculata (Pf) 

Lycium shockleyi (Ls) 

Yucca brevifolia (Yb) 

Salazaria mexicana (Sin ) 
Thamnosma montuna (Tin) 
llaplopappus cooperi (1 Ic) 


+ 0.34° 0 "Oh, + 0.29° Mp 
-0.28° La 

+ 0.75 Ys. + 0.29° Hr 
+ 0.27 °Lp, + 0.35° ° °Ta. + 0.63° 0 "Asp 
+ 0.20+ Kp, +0.28 "Hr 

+ 0.50"" "En. + 0.27 "Pfr, +0.25“Ca, + 0.29'Pf, +0.20" Ad 
-0.27°Ta. -0.25"Aea, +0.36“"°Pp 
+ 0.38°"Gs, +0.25"Oh, -0.00"°"Pp. -0.28° Aco 

+ 0.75"°"Ef 

+ 0.33'" "Mt. + 0.25 "Hs 
+ 0.50'" "Cs 

+ 0.33" "Ssp. +0.38'"La 
+ 0.08° "Sp 

+ 0.25'La. + 0.02' "Mp. +0.34'"As 
+ 0.33 "' Sa, +0.41" "Pfr, +0.75""IIs 
+ 0.27 Kp. + 0.41 0 °"Mt, +0.33 Hs, +0.03'"Ca, +0.82°°°Pf. 
+ 0.27°Tm 

+ 0.32°' "Ta, +0.25'Ssp 
+ 0.27 "Cl 

+ 0.29’Eh +0.28 "Ad 

+ 0.25"Sa. + 0.75°"Mt, +0.33 Pfr, +0.48"""Ta. +0.33'"'Ta 
+ 0.08'"Lt 

+ 0.33" "(is, + 0.25 "Cr 

+ 0.35""Cl.-0.27"Lt, +0.32 Cr, +0.18'"Ils, +0.01"°Asp 
+ 0.29° As, +0.02"'Oh 
-0.25'Lt 

+ 0.03" "Cl, + 0.33"'1 Is, + 0.01 * "Ta 
+ 0.3(7"'Lt 

+ 0.25’Kp, +0.03"'Pfr, +0.53”°Pf, +0.45°"Yb. +0.72° "Tin 
+ 0.29'Kp, + 0.82' "Pf, + 0.53' "Ca 
+ 0.35' "Sm 

+ 0.45°”Ca. +0.51" "Tin 
+ 0.35" °'Ls, +0.51'"'Tm 
+ 0.27'Pfr, + 0.72“ ° "Ca, +0.51°°°Yb 


’Denotes a correlalion al 5 percent probability level, ” at 1 percent probability, and ” * at 0,1 percent probability. 


1980 


Nevada Desert Ecology 


55 



Fig. 6. Special ordination on the first three axes derived from the principal component analysis. Species showing 
positive correlations are interconnected. Letters beside each circle indicate abbreviations for species' names. 
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Sp 
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MP 
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24. 

Ta 
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27. 
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= Artemisia spinescens 

11. Sa 
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28. 

Pp 

= Psorothamnus polydenius 

12. En 

= Ephedra nevadensis 

29. 

Ca 
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13. Gs 

= Crayia spinosa 

30. 

Pf 

= Primus fascicidata 

14. Lf 

= Lepidium fremontii 

31. 

Ls 

= Lycium shockleyi 

15. Oh 

= Oryzopsis hymenoides 

32. 

Sin 

= Salazaria mexicana 

16. Mt 

= Machaeranthera tortifolia 

33. 

Yb 

= Yucca brevifolia 

17. Pfr 

= Psorothamnus fremontii 

34. 

He 

= Haplopappus eooperi 
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points in an s-dimensional coordinate frame. 
Classification consists of subdividing the 
swarm of points into a number of disjointed 
sets. If the points chance to fall into several 
compact, widely separated groups, no diffi¬ 
culty arises and formal rules for effecting a 
classification are scarcely needed. According 
to Pielou (1969) this ideal result is rarely ob¬ 
tained when vegetation is randomly sampled. 
More often than not, the points representing 
the stands are diffusely scattered and any 
classification procedure is largely arbitrary. 

A way out of this difficulty is to ordinate 
the stands rather than to classify them. The 
purpose as in classification is still to simplify 
and condense the mass of raw data yielded by 
vegetation sampling in the hope that rela¬ 
tionships among the plant species, and be¬ 
tween them and the environmental variables, 
will be manifested. 

Ordination consists of plotting n stands in a 
space of fewer than s dimensions in such a 
way that none of the important features of 
the original s-dimensional pattern is lost. Or¬ 
dination has two great advantages over classi¬ 
fication. It obviates the necessity for setting 
up arbitrary criteria for defining the classes, 
and there is no need to assume that distinct 
classes, if there are any, are hierarchically re¬ 
lated. However, the compatability of the two 
approaches, viz., the classification and the or¬ 
dination techniques, was pointed out by An¬ 
derson (1965) and Goodall (1970), and in re¬ 
cent years it has become more common for 
classification and ordination to be used on 
the same data (Gray and Bunce 1972, Wil¬ 
liams and Walker 1974, Ayyad and El-Gho- 
nemy 1976). 

The importance of the study area from the 
phytogeographical point of view may be due 
to its position straddling the boundaries of 
the Great Basin to the north and the Mojave 
Desert to the south. The behavior of the biot¬ 
ic communities of the Mojave Desert in gen¬ 
eral, or in some of its sectors, has attracted 
the attention of many biologists. Shreve and 
Wiggins (1964) have described the Mojave 
Desert as showing its most distinctive devel¬ 
opment between 600 and 1200 m elevation 
(2000-4000 ft.). When it is followed thence 
toward the northeast or southeast, it loses 
some of its characteristic vegetational fea¬ 
tures and much of its distinctive flora. The 


basic structure of the vegetation throughout 
the Mojave Desert is open stands of L. triclen- 
tata and A. dumosa. On the western edge the 
plants are joined and to some extent replaced 
by A. tridentata, G. spinosa, T. axillaris, and 
other perennials; and at high elevations on 
the north, C. ramosissima, G. spinosa are 
dominant. 

In the northern sector of the Mojave 
Desert, particularly in the Nevada Test Site 
and its surroundings, phytosociological stud¬ 
ies have been carried out by many authors. 
Among these are Beatley (1963, 1969, 1974, 
1975, 1976), Allred et al. (1963), Rickard and 
Beatley (1965), Brown and Mason (1968), 
Wallace and Romney (1972), Romney et al. 
(1973), Tueller et al. (1974), and El-Ghonemy 
et al. (this volume). In these studies some of 
the vegetational units have been identified 
and named by various terms as communities, 
associations, types, subtypes, and vegetational 
groupings. In the most recent work by El- 
Ghonemy et al. (1980) the correlation be¬ 
tween the various vegetational groupings 
identified and the environmental variables 
have been demonstrated, as well as the vege¬ 
tation diversity and the successional trends 
among the different groupings. 

In the present study the application of the 
agglomerative clustering technique has 
proved useful in classifying stands into sever¬ 
al vegetational groupings. However, most of 
these groupings are not distinct. The mem¬ 
bers of each pair of groupings are, in most 
cases, linked together by having one or more 
of the dominant species in common. This, as 
Goodall (1954) mentioned, does not preclude 
the possibility of the classification for par¬ 
ticular ends, but it is generally more appro¬ 
priate to ordinate stands. 

In the present study the application of 
both Wisconsin and PCA ordination tech¬ 
niques emphasizes this idea and indicates that 
the vegetational groupings yielded by the 
clustering technique are generally inter¬ 
connected. However, the PCA technique has 
proved more efficient in segregating stands in 
a manner more or less similar to that 
achieved by the clustering approach. The 
better segregation of stands along the PCA 
ordination plan makes the description and 
correlation of the vegetational grouping 
more efficient. 
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The vegetational groupings identified in 
the study area are more or less similar to 
those previously identified by El-Ghonemy et 
al. (this volume) with stands clustered into 
groupings according to their leading domi¬ 
nant species (i.e., species with highest impor¬ 
tance values). Because in the multivariate 
analysis the similarity between stands in¬ 
volves the use of similarity functions that 
take into account the whole number of spe¬ 
cies involved in community structure, the 
stand’s characteristics for a given vegetation 
unit may not have the same leading dominant 
species. Consequently, the detailed structure 
of comparable vegetation groupings derived 
from the two approaches may not be the 
same. 

The spatial arrangement of stands along 
the different vegetational gradients provides 
evidence that variations in vegetation com¬ 
position are expressed by more than one axis 
of the ordination components. This implies 
that the distribution of vegetation in the 
study area is controlled by complexes of in¬ 
terrelated factors. These fall into three main 
groups. The first group relates to soil fertil¬ 
ity, as reflected by the concentrations of 
phosphorus, nitrogen, potassium, and other 
nutrient elements. The second complex ex¬ 
presses soil salinity, and the third complex re¬ 
lates to soil texture and water-retention ca¬ 
pacity. The role of soil fertility as a factor in 
the delimitation of the natural plant commu¬ 
nities has been stressed by several authors 
(Beadle 1954, 1962, El-Ghonemy 1966, 
Ayyad and El-Ghonemy 1976). 

In the present study the association of cer¬ 
tain vegetational groupings, for example, that 
of A. dumosa, with phosphorus-rich soil, has 
been demonstrated. The role of total soil 
salinity or of particular salt has been criti¬ 
cally reviewed by Chapman (1960), Gates et 
al. (1956), Ayyad and El-Ghareeb (1972), and 
Beatley (1976). The Atriplex communities at 
the Nevada Test Site have been described by 
Beatley (1976) as occupying sites at the high 
end of the soil salt content and fine particle 
gradients. These results are in bill agreement 
with the results obtained in the present study. 

Conclusions 

(1) The application of the multivariate 
analysis resulted in the segregation of the 


vegetational cover into more or less distinct 
groupings. 

(2) The application of the principal com¬ 
ponent analysis for the ordination of the veg¬ 
etation data resulted in the segregation of 
stands into more or less distinct sets, com¬ 
parable in structure to the vegetation group¬ 
ings derived from the cluster analysis. 

(3) The physical and chemical properties 
of the soil play a definite role in the delinea¬ 
tion of the vegetational groupings. 

(4) The ordination of species scores may 
result in groupings of correlated species that 
could be floristically different from those ob¬ 
tained if ordination is based on stand scores. 

Acknowledgments 

The multivariate analyses were carried out 
using a Burroughs B-6700 at the Utah State 
University Computer Center. The authors are 
indebted to Dr. W. Valentine for his assis¬ 
tance in data processing. 

Support for this study was hirnished in 
part by Contract EY-76-C-03-0012 between 
the U.S. Department of Energy and the Uni¬ 
versity of California and an IBP Desert 
Biome subcontract with Utah State Univer¬ 
sity and the University of California. 

Literature Cited 

Allred, D. M„ D. E. Beck, and C. D. Jorgenson. 1963. 
Biotic communities of the Nevada Test Site. Brig¬ 
ham Young Univ. Sci. Bull., Biol. Ser. 2(2):52. 
Anderson, D. J. 1965. Classification and ordination in 
vegetation science: controversy over a nonex¬ 
istent problem? J. Ecol. 53:521-526. 

Ayyad, M. A., and R. El-Ghareeb. 1972. Micro¬ 
variations in edaphic factors and species distribu¬ 
tion in a Mediterranean salt desert. Oikos 
23(1): 125-131. 

Ayyad, Nl. A., and A. El-Ghonemy. 1976. Phytosociolo- 
gical and environmental gradients in a sector of 
the western desert of Egypt. Vegetation 
31(2):93-102. 

Beadle, N. G. W. 1954. Soil phosphate and delimitation 
of plant communities in eastern Australia. 1. 
Ecology 35:370-375. 

_ 1962. Soil phosphate and delimitation of plant 

communities in eastern Australia. II. Ecology 
43:281-344. 

Beals, E. W. 1960. Forest bird communities in the 
Apostle Islands of Wisconsin. Wilson Bull. 
72:156-181. 

Beatley, J. C. 1963. Vegetation and environment of the 
Nevada Test Site. Ecol. Soc. Amer., Bull, (abstr.) 
44:123. 


58 


Great Basin Naturalist Memoirs 


No. 4 


_ 1969. Vascular plants of the Nevada Test Site, 

Nellis Air Force Range, and Ash Meadow (north¬ 
ern Mojave and southern Great Basin deserts, 
south central Nevada). USAEC Report, UCLA 
12-705. 

. 1974. Effect of rainfall and temperature on the 
distribution and behavior of Larrea tridentata 
(creosote bush) in the Mojave Desert of Nevada. 
Ecology 55:245-261. 

. 1975. Climate and vegetation pattern across the 
Mojave Great Basin Desert transition of southern 
Nevada. Amer. Midi. Natur. 93(l):53-70. 

_ 1976. Vascular plants of the Nevada Test Site and 

central southern Nevada. ERDA Technical Infor¬ 
mation Center. T1D-26851. 803 p. 

Bray, J. R., and J. T. Curtis. 1957. An ordination of the 
upland forest communities of southern Wiscon¬ 
sin. Ecol. Monogr. 27:325-349. 

Brown, K. W., and B. J. Mason. 1968. Range survey 
area 18, Nevada Test Site. U.S. Dept, of Health, 
Education and Welfare, Public Health Service. 
SWR116-52. 

Chapman, V. J. 1960. Salt marshes and salt deserts of 
the world. Interscience Pub., Inc., New York. 352 

P- 

El-Ghonemy, A. A. 1966. Soil-vegetation relationships 
for some major plant communities in south¬ 
western New South Wales, Australia. Unpub¬ 
lished dissertation. University of New England, 
Armidale. 

El-Ghonemy, A. A., A. Wallace, and E. M. Romney. 
1980. Socioecological and soil-plant studies of the 
natural vegetation in the northern Mojave 
Desert-Great Basin desert interface. Great Basin 
Nat. Mem. 4:71-86. 

Gates, D. H., A. L. Stoddart, and W. C. Cook. 1956. 
Soil as a factor influencing the plant distribution 
on salt deserts of Utah. Ecol. Monogr. 
26:155-173. 


Goodall, D. W. 1954. Objective methods for the classi¬ 
fication of vegetation. III. An essay on the use of 
factor analysis. Aust. J. Bot. 2:304-324. 

. 1970. Statistical plant ecology. Ann. Rev. Ecol. 
Syst. 1:99-124. 

Gray, E. J., and G. H. Bunce. 1972. The ecology of 
Morecambe Bay. VI. Soils and vegetation of salt 
marshes: A multivariate approach. J. Appl. Bot. 
9:111-234. 

Pielou, E. C. 1969. An introduction to mathematical 
ecology. Wiley (Interscience), New York. 

Rickard, W. H., and J. C. Beatley. 1965. Canopy cov¬ 
erage of the desert shrub vegetation mosaic of 
Nevada Test Site. Ecology 47:524-529. 

Romney, E. M., V. Q. Hale, A. Wallace, O. R. Lunt, 
J. D. Childress, H. Kaaz, G. V. Alexander, J. E. 
Kinnear, and T. L. Ackerman. 1973. Some char¬ 
acteristics of soil and perennial vegetation in 
northern Mojave Desert areas of the Nevada Test 
Site. AEC Report UCLA 12-915, Biomed. and 
Environ. Res. Report TID-4500. 

Shreve, F., and W. Wiggins. 1964. Vegetation and flora 
of the Sonoran Desert. Vols. 1 and 2. Stanford 
University Press, Stanford, Calif. 

Sneath, P. H. A., and R. R. Sokal. 1973. Numerical tax¬ 
onomy. Freeman, San Francisco. 

Tueller, P. T., D. B. Allen, R. Everett, and J. B. 
Davis. 1974. The ecology of Hot Creek Valley, 
Nevada, and nonradiation effects of an under¬ 
ground nuclear detonation. U.S. Atomic Energy 
Commission, R 96, NVO 409-2. 

Wallace, A., and E. M. Romney. 1972. Radioecology 
and ecophvsiologv of desert plants at the Nevada 
Test Site. National Technical Information Ser¬ 
vices, USAEC Report T1D-25954. 

Williams, C. H., and H. B. Walker. 1974. The vegeta¬ 
tion of tropical African Lake: classification and 
ordination of the vegetation of Lake Chilwa (Ma¬ 
lawi). J. Ecol. 62(3):831-853. 




